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Abstract — This paper addresses the synchronization rate of 
weakly connected pulse-coupled oscillators (PCOs). We prove 
that besides coupling strength, the phase response function is also 
a determinant of synchronization rate. Inspired by the result, we 
propose to increase the synchronization rate of PCOs by design- 
ing the phase response function. This has important significance 
in PCO-based clock synchronization of wireless networks. By 
designing the phase response function, synchronization rate is in- 
creased even under a fixed transmission power. Given that energy 
consumption in synchronization is determined by the product of 
synchronization time and transformation power, the new strategy 
reduces energy consumption in clock synchronization. QualNet 
experiments confirm the theoretical results. 

Index Terms — Synchronization rate, pulse-coupled oscillators, 
phase response function, wireless networks 

I. Introduction 

In recent years, synchronization of oscillating dynamical 
systems is receiving increased attention. One particular class 
of oscillating dynamical systems, pulse-coupled oscillators 
(PCOs), are of considerable interest. 'Pulse-coupled' means 
that oscillators interact with each other using pulse-based 
communication, i.e., they can achieve synchronization via the 
exchange of simple identical pulses. PCO has been used to 
describe many biological synchronization phenomena such 
as the flashing of fireflies, the contraction of cardiac cells, 
and the firing of neurons [1]. Recently, with the progress 
of ultra-wide bandwidth (UWB) impulse radio technology 
[2], the PCO based synchronization scheme has also been 
applied to synchronize wireless networks [3], [4], [5], [6], [7]. 
Since it is implemented at the physical layer or MAC (Media 
Access Control) layer, it eliminates the high-layer intervention. 
Moreover, message exchanging in PCO-based synchronization 
strategy is independent of the origin of the pulses, which 
avoids requiring memory to store time information of other 
nodes [3]. Therefore, the PCO-based synchronization scheme 
has received increased attention in the communication com- 
munity recently. 
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Despite considerable work on synchronization conditions, 
there remains a lack of research on the synchronization rate for 
PCOs, especially for PCOs with a general coupling structure 
other than the commonly studied all-to-all structure. The 
synchronization rate is crucial in synchronization processes 
[8]. For example, in the clock synchronization of wireless 
networks, the synchronization rate is a determinant of the 
consumption in energy, which is a precious system resource 
[3]. 

This paper analyzes the synchronization rate of weakly 
connected PCOs in the presence of combined global cues (also 
called leader, or pinner in the language of pinning control 
[9]) and local cues (alternatively, local coupling). The network 
structure is considered because in the clock synchronization of 
wireless networks, usually different time references are syn- 
chronized through internal interplay between different nodes 
and external coordination from a global time base such as 
GPS [10]. Due to pulsatile coupling, the synchronization rate 
of PCOs are very difficult to study analytically [1]. Based on 
the assumption of 'weakly connected', we study the problem 
using phase response functions. A phase response function 
describes the phase correction of an oscillator induced by a 
pulse from neighboring oscillators or external stimuli [11]. 
Under the assumption of weak coupling, we transform the 
PCO model into a simpler phase model, based on which we 
analyze the synchronization rate of PCOs. In fact, we will 
prove that the synchronization rate is determined not only by 
the strength of global and local cues, but also by the phase 
response function. This means that different phase response 
functions bring different synchronization rates even when the 
coupling strength is fixed. This has great significance for 
synchronization strategies such as the clock synchronization 
of wireless networks, where the phase response function is a 
design parameter. By designing the phase response function, 
we increase the synchronization rate under a fixed coupling 
strength. Given that the total energy consumption in synchro- 
nization is determined by the product of synchronization time 
and transmission power (corresponding to coupling strength 
and topology) [3], [12], the new strategy reduces energy 
consumption in synchronization in that synchronization time 
is reduced under a fixed transmission power. It is worth 
noting that the assumption of weak coupling is well justified 
by biological observations: the amplitudes of postsynaptic 
potentials are around 0.1 mV, which is small compared with 
the amplitude of excitatory postsynaptic potential necessary to 
discharge a quiescent cell (around 20 mV) [13]. In PCO-based 
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wireless network synchronization schemes, weak coupling is 
also necessary to guarantee a robust synchronization [3]. 

II. Problem formulation and Model 

TRANSFORMATIONS 

Consider a network of N pulse-coupled oscillators, which 
will henceforth be referred to as 'nodes'. All oscillator nodes 
or a portion of them can receive alignment/entrainment in- 
formation from an external global cue (also called leader, or 
pinner in the language of pinning control [9]). 

We denote the dynamics of the oscillator network as 

%g = fg( x g) 

x i =fi(x i ) + g i S(t-t g )+l (1) 

l<j<N,j& 

for i = 1,2,..., JV", where x g G [0, 1] and xi G [0, 1] denote 
the states of the global cue and oscillator nodes, respectively. 
f g and /, describe their dynamics, g^ > denotes the effect of 
the global cue's firing on oscillators v. when x g reaches 1 (at 
time instant t g ), it fires and returns to 0, and at the same time 
increases oscillator i by an amount g^. I > and G {0, 1} 
denote the effect of oscillator j's firing on oscillator i: when 
Xj reaches 1 (at time instant tj), it fires and resets to 0, and 
at the same time pulls oscillator i up by an amount lcti ; j. The 
increased amount is produced by dirac function S(t), which is 
zero for all t except t = and satisfies /_ 5(t)dt = 1. 

Remark 1: If (or a, ,) is 0, then oscillator i is not affected 
by the global cue (or oscillator j). 

Assumption 1: We assume a^j = aji, which is common in 
wireless networks [14], [15]. 

Assumption 2: We assume weak coupling [13], i.e., g and 
I satisfy j«l and ! < 1. 

Assumption [2] follows from the fact that the amplitudes of 
postsynaptic potentials measured in the soma of neurons are 
far below the amplitude of the mean excitatory postsynaptic 
potential necessary to discharge a quiescent cell [13], it is 
also required in PCO-based wireless network synchronization 
strategies to ensure the robustness of synchronization [3]. 

Based on Assumption 12 the system in (H) can be described 
by the following phase model using the classical phase reduc- 
tion technique and phase averaging technique [11], [16]: 




l<j<N,j& 

(2) 

for i = 1,2,..., N, where g G [0, 2tt) and 9 t G [0, 2tt) 
denote the phases of the global cue and oscillator i, respec- 
tively. Q g (x) and Qi(x) are phase response functions and are 
often referred to as phase response curves in biological study. 
They are periodic functions with period 2tt [11], [16]. T is 
the period of the global cue. w g and Wi denote the natural 
frequencies of the global cue and oscillator i, respectively. 

Remark 2: The transformation from (Q]i to (O is a standard 
practice in the study of weakly connected PCOs and it is 
applicable to any limit-cycle oscillation function fa and f g 



[11]. The detailed procedure has been well documented in 
[16], Chapter 9 of [13], and Chapter 10 of [11]. 

Assumption 3: In the paper, we assume that Q p (p = 
{I, <?}) satisfy the following conditions: 

Q p (0) - 0; Vz G (-7T, tt), > 0; Q p (-x) = -Q P (x) 

(3) 

Remark 3: Assumption [3] gives an advance-delay phase 
response function, which is common in biological oscillators 
[11]. Moreover, given that in wireless networks, the phase 
response function is a design parameter, Assumption [3] will 
simplify analysis and design, and as shown later, such phase 
response functions will also lead to good synchronization 
properties. 

Solving the first equation in (O gives the dynamics of the 
global cue 9 g = w g t + £ s , where the constant £ g denotes the 
initial phase of the global cue. To study if local oscillators can 
be synchronized to the global cue, it is convenient to study the 
phase deviation of local oscillators from the global cue. So we 
introduce the following change of variables: 

0i = 6g + £ = Wgt + ^+ & (4) 

Therefore G [— tt, tt] denotes the phase deviation of the ith 
oscillator from the global cue. Substituting into (fJI yields 
the dynamics of phase deviations 

ii = A i -^Q g (Z i ) + i a idQi& - ti) ( 5 ) 

for i = 1, 2, . . . , N, where A, = w, — w g . In (0, the oddness 
property of function Q g is exploited. 

Assumption 4: In this paper, we assume Wi = w g is 
satisfied for all i — 1,2, . . . , N, i.e., all the oscillators have 
the same natural frequency as the global cue. 

Using Assumption @] (0 reduces to: 

& = -§Q s &) + ^ £ a id Q t (Zi-Si) (6) 

for£=l,2,...,JV. 

Thus far, by analyzing the properties of ©, we can obtain 
the roles of global and local cues as well as phase response 
functions in the synchronization of PCO networks: 

• Synchronization: If all £j asymptotically converge to 0, 
then we have 9\ = 82 = ■ ■ ■ = On when time goes to 
infinity, meaning that all the nodes are synchronized to 
the global cue. 

• Exponential bound on the synchronization rate: From 
dynamic systems theory [17], the synchronization rate is 
determined by the rate at which £j decays to 0, namely, 
it can be measured by the maximal value of a (a > 0) 
satisfying the following inequality for some constant C: 

U(t)\\<Ce- at U(0)\\, f=[6 £2 ... 6v ] T 

(7) 

where || • || is the Euclidean norm. A larger a leads to a 

faster synchronization rate. 
Assigning arbitrary orientation to each interaction, we can 
get the N X M incidence matrix B (M is the number of non- 
zero ciij (1 < i < N, j < i), i.e., the number of interaction 



3 



edges) of the interaction [18]: 5jj — 1 if edge j enters node Note that (-B T £)i (1 < i < M) are in the form of 



i, B t 



1 if edge j leaves node i, and Bij = otherwise. 



Then using graph theory, we can write (O in a more compact 
matrix form: 



(8) 



where £ is given in ©, G = diagjgi, g 2 , . .., ffjv}, and 
diag{»} denotes a diagonal matrix with elements {•} on the 
diagonal. 

III. Synchronization of Pulse Coupled 
Oscillators 



" 2 ' 2 > 



A. When all £j are within (- 

Theorem 1: For the oscillator network in ([8j, if all £, are 
within [— e, e] for some e G [0, J), then the network synchro- 
nizes to the global cue when at least one gi is positive and the 
local coupling topology a^j is connected. Here 'connected' 
means that there is a multi-hop path (i.e., a sequence with 
nonzero values a ijTOl , a mi , TO2 , a mp _ ump , a mp ,j) from 
each node i to every other node j. 

Proof: We first prove that for any < e < when 
£ G [— e, e] x ... x [— e, e] = [— e, e]^ where x is Cartesian 
product, they will remain in the interval, i.e, [— e, e] is 
positively invariant for ([8]). To this end, we only need to check 
the direction of the vector field on the boundaries. If = e, 
we have < — £j < 2s < ir for 1 < j < N, so from © 
and the properties of phase response functions in Assumption 
[3] fj < holds. Hence the vector field is pointing inward in 
the set, and no trajectories can escape to values larger than 
s. Similarly, we can prove that when £j = — e, £j > holds. 
Thus the vector field is pointing inward in the set, and no 
trajectories can escape to values smaller than — e. Therefore 
[— s, s] is positively invariant for any < e < ^. 

Next we proceed to prove synchronization. Construct a 
Lyapunov function as V = ^£ T £- V is non-negative and will 
be zero if and only if all are zero, meaning that all oscillators 
are synchronized to the global cue. 

Differentiating V along the trajectories of dHJ yields 



V 



ft 



±eGs 1 z-! t eBs a B T z 
i 



T 



f (GSi + IBS 2 B T ) £ 



where Si and S 2 are given by 



diag 



4 



N 



(9) 



(10) 



(S T £)i 



Af 



(11) 



with (B T £)i (1 < i < M) denoting the ith element of the 
M x 1 dimensional vector B T £. 

According to dynamic systems theory [17], if GSi + 
IBS 2 B T in (0 is positive definite, then V is always negative 
when £ 7^ and V will decay to zero exponentially, meaning 
that £ will converge to zero and all oscillators are synchronized 
to the global cue. 



Cm ~ €n (1 < m,n < N), it follows that (B T £)i are 
restricted to [— 2e, 2e] when all £j are in [— e, e] for some 
£ G [0, 77). Given that in (— 71", 7r), Q g (x) and Qi(x) satisfy 



>0, 



> 0, it follows that 5i and 5 2 are positive 



definite, and thus the following inequalities are satisfied for 
some positive constants o\ and a 2 : 

Si > ail, S 2 > a 2 I, 

Qg{x) . Ql(x) 7T 

o\ = mm — - , g 2 = mm , e 6 0, — ) 

-s<x<e X -2e<x<2e X 2 

(12) 

So we have GSi + IBS 2 B T > a x G + a 2 lBB T , which in 
combination with (O produces 
1 



V<--e {a.G + ^lBB 1 )^ 



(13) 



Next we prove that o\G + a 2 lBB T is positive definite, which 
leads to V < for £ ^ 0. 

It can be easily verified that a\G + a 2 lBB T is of the 
following form: 

o-i G + a 2 lBB T = <Tidiag{3i, g 2 , gw} + o 2 lL (14) 

with L E jiNxn cons tructed as follows: for i ^ j, its 
(i,j)th element is —Qij, for i = j, its (i,j)th element is 

JV 

ai.m- Since o\, a 2 , and I are positive, and gi, are 

rn — l,m^i 

non-negative, it follows from the Gershgorin Circle Theorem 
that <j\G + a 2 lBB T only has non-negative eigenvalues [19]. 
Next we prove its positive definiteness by excluding as an 
eigenvalue. 

Since the topology of local coupling a^j is connected, 
o\G + a 2 \BB T is irreducible according to graph theory [19]. 
This in combination with the assumption of at least one 
non-zero guarantees that o^G + a 2 lBB T is irreducibly 
diagonally dominant. So from Corollary 6.2.27 of [19], we 
know the determinant of a\G+a 2 lBB T is non-zero and hence 
is not its eigenvalue. Therefore o\G + a 2 lBB T is positive 
definite, and V will converge to 0. ■ 

B. When all are within (— 7r, 71") and the maximal/minimal 
is outside (— f , f) 

Theorem 2: For the oscillator network in (0, if all £j are 
within [—£,£] for some e G [f , tt) and the maximal/minimal 
£i is outside (—5, ■§), then the oscillator network will syn- 
chronize to the global cue when all nodes are connected to 
the global cue, and the following relations are satisfied: 

CT4?A max (i?-B T ) 



> 



03 



9i>— ^ a i,jl2, i = 1,2 N 

7i 



(15) 



l<3<N,j& 

where A max denotes the maximal eigenvalue, g n 
min{5i, g 2 , g N }, and 

Qi(x) 



Qg{x) 

(73 = mm — , 04 = max 

-£<x<s X -2e<x<2e 



7x = mm Q g {x), j 2 = max -Qi{x), e G [-, tt) 

0<x<e 0<x<2s Z 



(16) 
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Proof: Following the line of reasoning of Theorem Q] we 
can prove that if the second inequality in ( Tl3T > holds, then for 
any ^ < e < tt, [—e, e] is positively invariant for (O, i.e., 
for £ G [— e, e] N , it will always remain in the interval. Next 
we proceed to prove synchronization. 

Choose the same Lyapunov function V as the proof of 
Theorem Q] Then we have 

V = i T i = -^fGSit - U T BS 2 B T £, 

\ 1 (17) 

= --£ T (GSi + IBS 2 B T ) i 

where Si and 52 are given in (TTOb and (fTTT i. respectively. 

When the maximal/minimal £j is outside (— §), Cm ~~ 
C« (1 < m, n < N) may be outside (— tt, tt). So in (fTOb 
and (fTTT i. the domain of Q g {x) is within (— 7r, 7r), on which 
Q g {x) satisfies Q g (x)/x > 0, and the domain of Qi(x) is 
not restricted to (—tt, tt), outside of which, Qi(x)/x may be 
positive or negative. Therefore Si is still positive definite, but 
£2 may be positive definite, negative definite or indefinite. 
From the definition of 03 and 04 in ( TTSI l. we have: 

GSi > a 39min , -IBS 2 B T < a i l\ max (BB T ) 

Notice that Qi(x) is periodic with period 2tt, and 2li£i > 
holds for all — tt < x < tt, we know for any | < e < tt, if 
x G [— 2e, — 7r], then Q;(a;o) = Qz(27r + a:o) > holds since 
2tt + Xq resides in the interval [2(tt — e), tt]. Thus it follows 

QjXxoX = Q^Xo+gTT) < Q which means > Q 

Therefore, (113b guarantees the positive definiteness of 
GSi + IBS 2 B T , and hence the synchronization of the os- 
cillators to the global cue. ■ 

Remark 4: Theorem [2] indicates that when the maxi- 
mal/minimal phase difference is outside (—5, ■§), all oscilla- 
tors have to connect to the global cue to ensure synchroniza- 
tion to the global cue. This is consistent with existing results 
which have shown that for some initial conditions (even with 
measure zero), PCOs cannot be synchronized by local coupling 
[1]. In fact, most of the existing results on PCOs are based on 
all-to-all connection, which amounts to g m j n > 0. 

Remark 5: Theorem [2] reveals that a strong local cue does 
not necessarily benefit synchronization when phase difference 
is outside (— ■§, § ). This is also consistent with [20] which 
shows that synchronization may not be achieved despite arbi- 
trarily strong local coupling. 

IV. Synchronization Rate of Pulse Coupled 
Oscillators 

Based on a similar derivation, we can get a bound on the 
exponential synchronization rate: 

Theorem 3: For the oscillator network in dH), define ai, a 2 
as in (fl2l i, and 03, 04 as in ( fl6] l, then 

• when all £j are within [— e, e] for some e G [0, and the 
conditions in Theorem 1 hold, the synchronization rate is 
no worse than 

ai = min {£ T (01G + a 2 lBB T ) £/(£ T 0} /T 

f (18) 
= A min (<7iG + a 2 WB T ) IT 



• when the maximal/minimal & is outside (— f , f ) and the 
conditions in Theorem 2 hold, the synchronization rate is 
no worse than 

Oi 2 = (<J3g m in - C4^A 

max 

(BB T )) IT (19) 
Proof: First consider the case that all £j are within [— e, e] 
for some e G [0, |). From (O, it follows 

V" < -2ai^ = -2c*iV (20) 
with o?! defined in ( fT8l ), which further means that 

V(t) < C 2 e~ 2ait V(0) U(t)\\ < Ce- ait U{0)\\ (21) 

holds for some positive constant C. Thus the synchronization 
rate is no worse than cti in ( fTST l. 

Similarly, we can prove from ( fTTI i that when the maxi- 
mal/minimal & is outside (— § ), the following relationship 

V(t) < C 2 e- 2Q2 V(0) \\Z(t)\\ < C7e^ Q2t ||e(0)|| (22) 

holds for some positive constant C. Thus the synchronization 
rate is no worse than a 2 in (1191 . ■ 

Remark 6: When Q g and are sinusoidal functions, all 
£i are constrained in the interval (— f , and there is no 
global cue (G = 0), using 9 = X) W as re f erence > we 

l<i<N 

can define £j as ^ = ^ — 8. Since £ T 1 = holds for 
1 = [1, 1, . . . , 1] T , the constraint £ T 1 = is added to the 
optimization mm {£ T (aiG + a 2 lBB T ) C/(^ T 0} l T in ©■ 

Given that G = and BB T is the Laplacian matrix of 
interaction graph and hence has eigenvector 1 with associated 
eigenvalue [19], A m ;„ in ( fT8l reduces to the second smallest 
eigenvalue, which is the same as the convergence rate in 
section IV of [21] obtained using contraction analysis. 

From application point of view, it is important to analyze 
how synchronization rate is affected by the phase response 
function and the strengths of global and local cues. According 
to ( fl6l l and ( fl9l l. it is clear that the synchronization rate 
increases with an increase in g m i n and ^ g ■ But how the 
phase response function and the strength of the global cue 
affect the synchronization rate when all £j are within (— f , f) 
is not clear. (In this case, <? m i n may be zero since some 
oscillators may not be connected to the global cue.) In fact, 
we can prove that in this case the synchronization rate also 

O (x) 

increases with an increase in ; and the strength of the 
global cue: 

Theorem 4: The synchronization rate of ([8]) increases with 
an increase in the strength of the global cue. It also increases 
with an increase in ^ g . 

X 

Proof: As analyzed in the paragraph above Theorem 
[4] we only need to prove the statement when all ^ G 
[— e, e] for some e G [0, ^), i.e., eei in ( TT8l is an in- 
creasing function of and SaljEl Recall from dT4b that 
01G + a 2 \BB T is an irreducible matrix with non-positive 
off-diagonal elements, so there exists a positive fj, such that 
/.il — (aiG + a 2 lBB T ) is an irreducible non-negative ma- 
trix. Therefore, A max [pi — [piG + a 2 lBB T )) is the Perron- 
Frobenius eigenvalue of [il — (criG + a 2 lBB T ) and is positive 
[19]. Given that for any 1 < i < N, /j, - A^ctiG + a 2 lBB T ) 
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is an eigenvalue of matrix fil — {a±G + a 2 lBB T ) where \ 
denotes the ith eigenvalue, we have 

Am»(/i/ -(a x G + a 2 lBB T )) 



fi-\ m in(viG + <r 2 lBB T ) 



i.e., 



0L\ 



\ min (<JiG + a 2 lBB T ) 

a l G + G 2 lBB T )) 



/' 



T 

(/// 



T 



Given that the largest eigenvalue (also called the Perron- 
Frobenius eigenvalue) of fil — (a\G + a 2 lBB T ) is an in- 
creasing function of any of its diagonal element [19], which 
is a decreasing function of gi and Sai£2 ; it follows that 
Amax (fil — (ciG + a 2 WB T )^j is a decreasing function of 
both gi and < ^ a( - x - ) , meaning that ot\ is an increasing function 
of 9l and 2^1. ■ 

Remark 7: The role of the local cue is not discussed in 
Theorem [4] In fact, the role of the local cue depends on 
the value of when all are within (— f , §), S 2 in 
dTTT > is positive definite, so ^ T BS 2 B T ^ in (O is positive, 
meaning that the local cue will increase the synchronization 
rate. Whereas when the maximal/minimal £j is outside of the 
interval (— ?, S 2 in <HD can be positive semi-definite, 
negative semi-definite or indefinite, £ T BS 2 B T !; in ( fTTI i can 
be positive, negative or zero, thus the local cue may increase, 
decrease or have no influence on the synchronization rate. This 
conclusion is confirmed by QualNet experiments in Sec. |VI] 

Remark 8: From Theorem [3] and Theorem HI one can see 
that in addition to the strength of coupling, i.e., gi and 
I, the phase response function Q g also influences the syn- 
chronization rate. This has significant ramifications for the 
clock synchronization of wireless networks using PCO-based 
strategies [3], [6], [12], where the phase response function is 
a design parameter: the synchronization rate can be increased 
by choosing appropriate phase response functions, even with 
transmission power (corresponding to coupling strength and 
topology) fixed, therefore leading to a reduced energy con- 
sumption. This will be addressed in Sec. [V] 

V. Design of Phase Response Functions 

As stated in Sec. IIVI the phase response function is an 
important determinant of the synchronization rate of PCOs. 
This has important ramifications for PCO-based synchroniza- 
tion strategies of wireless networks, where the phase response 
function is a design parameter. 

PCO-based synchronization strategies are attracting in- 
creased attention in the communications literature [3], [6], [7], 
[12]. As with most synchronization strategies, a network using 
a PCO-based synchronization strategy makes a distinction 
between an acquisition stage where the network synchro- 
nizes and the communication stage where nodes transmit and 
receive data [14]. In PCO-based synchronization strategies, 
every node of the network acts as a PCO, nodes interact 
through transmitting replicas of a pulse signal, which can be 
a monocycle pulse in a UWB network [3] or preambles in 
IEEE 802.11 networks [22]. PCO-based strategies have many 



advantages over conventional synchronization strategies [3]: 
they are implemented at the physical layer or MAC layer, 
which eliminates the high-layer intervention; their message 
exchanging is independent of the origin of the signals, which 
avoids requiring memory to store time information of other 
nodes. 

In all existing PCO-based synchronization strategies, the 
oscillator model is directly adopted from a biological source, 
leading to a fixed phase response function. Our finding sug- 
gests that even under a fixed coupling strength, one can 
increase the synchronization rate by designing the phase 
response function. This can reduce energy consumption in 
clock synchronization since the total energy consumption in 
a synchronization process is determined by the product of 
transmission power (corresponding to coupling strength and 
topology) and the time to synchronization. Next, we show 
that by designing the phase response function, we can indeed 
increase the synchronization rate. 

We focus on a class of phase response functions in tanh 
form, for reasons outlined below: 



Q(x) 



tanh(x/e) x 



when x € [— tt, tt] (23) 



tanh(7r/e) 

where e > is a free parameter and will be designed to 
achieve a faster synchronization rate. 

Fig. Q] gives the plot of Q(x) in ( l23l . Since Q(x) is 2ir- 
periodic, only its value in the interval [0, 2ir] is plotted. Q(x) 
is an advance-delay phase response function, i.e., external 
pulsatile input either delays or advances an oscillator's phase, 
depending upon the timing of input. Advance-delay phase 
response functions have been widely used in the biology 
community: they can well characterize the dependence of 
neurons' response to small depolarizations, i.e., an excitatory 
postsynaptic potential (EPSP) received after the refractory 
period delays the firing of the next spike, while an EPSP 
received at a later time advances the firing. The most widely- 
used neuron model, i.e., the Hodgkin-Huxley model, also has 
advance-delay phase response functions [23]. From Theorem 




Fig. 1, The shape of tanh type advance-delay phase response functions in 
equation i23\ . 



21 we know in addition to the strength of the global and 
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local cues, the phase response function of the global cue also 
determines the synchronization rate: the larger Salg) i Sj me 
faster the synchronization rate. In the following, we will show 
that the synchronization rate can be increased by designing e, 
a parameter in the phase response function. 

Theorem 5: For the pulse-coupled oscillator network in ([H) 
and the phase response function of the global cue in the 
form of d23l ). if all £j are within (— tt, tt), the synchronization 
rate increases with a decrease in e, no matter whether the 
maximal/minimal is within or outside (— f-, §). 

Proof: According to Theorem|H the synchronization rate 
increases with an increase in 2sl£2. So we only need to prove 
that fizlsl increases with a decrease in e for x g (— 7T, tt). 

Using the tanh type phase response function, we have 

Q g (x) tanh(x/e) 1 



x *tanh (7 r/ £ ) «> * € ^ (24) 

Since j n (F24I is a smooth function of e, we can calculate 

its derivative with respect to e: 



d(^£M) 7rsech z (^)tanh(^) - x sech^) tanh(^) 



,2/7T 



cfe 



xe 2 tanh 2 (^) 



(25) 



To prove that a ' increases with a decrease in e, we need to 



prove that — rf g — - is negative. Since e tanh ( j) is positive, 
we only need to prove that d26l ) is negative for —tt < x < 7r: 

— (tt sech 2 (— ) tanh( — ) — x sech 2 (— ) tanh( — )] (26) 

a; V e e e e / 

Using properties of hyperbolic functions, we can rewrite 
as follows: 







it e e 



(27) 



hence the problem is reduced to proving the negativity of 
f(x, e) in d28| i for — tt < x < tt. 



f(x,e) = 



tt (e~* — e e — x (e 2 ? — e * J 



.r 



(28) 



Using a Taylor expansion, equation (128t can be further 
rewritten as 



/(x,e) = - 
x 



, 2x 
2tt — 



(f) 3 , ) 5 
3! 



5! 



-2a; — 



3! 5! 

2 / Trg r-x(^r ^e-ff-xc-ff 

x\ 3! 5! 



(29) 



which is negative for all x in — 7r < x < tt. 



So 



is negative, thus 



increases with a decrease 



in e, which completes the proof. 



VI. QualNet Experiments 

We use a high-fidelity network evaluation tool (QualNet) 
to illustrate the proposed strategy. QualNet is a commercial 
network platform that has been widely used to predict the 
performance of MANETs, satellite networks and sensor net- 
works, among others [24]. 

In the implementation, we constructed a wireless network 
composed of 19 nodes (including 1 global cue). Each node 
has a counter as clock and stores the phase response function 
(as shown in Fig. [TJ in a lookup table. Upon receiving a 
pulse, a node shifts its phase by an amount determined by 
its current time and the phase response function in the lookup 
table. The structure of the network is illustrated in Fig. |2] 
where node number is the global cue. A broadcasting-based 
MAC layer protocol is adopted to establish the pulse based 
communication between different nodes, which is represented 
by the circles in Fig. |2] Although a broadcasting scheme is 
used, the communication in the network is not all-to-all due 
to limited transmission range. The interaction topology is also 
illustrated in Fig. [2] which can be verified to be a connected 
graph. 
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Fig. 2. Schematic of the network (left) and interaction matrix (right) used in 
the QualNet implementation. In the interaction matrix, ' 1 ' denotes interaction 
between two nodes, and '0' denotes no interaction. 



We first considered the case where £j G (—§,?)■ From 
Theorem Q] all nodes can synchronize to the global cue if 
at least one gi is non-zero. To confirm the prediction, we 
only connected oscillator 1 to the global cue with g\ = 0.01. 
We set the natural frequency to 2tt, the strength of local 
coupling to I = 0.01, and implemented the PCO network under 
different phase response functions in (l23l l. i.e., different pairs 
of e g (e in Qg) and ei {e in Qi). The network was synchro- 
nized, confirming TheoremQ] To illustrate our phase-response- 
function based design strategy, we recorded the average time 
to synchronization from 100 runs for each pair of e g and 
In each run, we set the initial value of the global cue as 
and chose the initial values of local nodes randomly from a 
uniform distribution on (— J). The times to synchronization 
are given by the first element of each 2-tuple in Table I. We 
can see that the synchronization rate increases with a decrease 
in e g , which confirms the theoretical results in Sec. [V] The 
total energy consumption in the synchronization process is 
also recorded and averaged over the 100 runs. The results are 
given by the second element of each 2-tuple in Table I. With 
a decrease in e g , the energy consumption indeed decreases, 
which confirms the effectiveness of our design methodology 
in Sec. El 
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TABLE I 

Time to synchronization [s](1st element of each 2-tuple) and energy consumption [10 _3 J](2nd element of each 2-tuple) under 

DIFFERENT PHASE RESPONSE FUNCTIONS (fj g (- £ , £ ), I = gi = 0.01, 92 = 93 = • ■ • = 9JV = 0) 



e g \ £l 


0.05 


0.1 


0.2 


0.4 


0.8 


1.6 


0.4 


(37.23, 729.39) 


(36.76, 719.99) 


(36.55, 715.79) 


(36.23, 709.05) 


(35.83, 701.39) 


(43.37, 852.19) 


0.8 


(37.88, 742.39) 


(37.40, 732.79) 


(37.88, 742.39) 


(36.73, 719.39) 


(36.44, 713.59) 


(44.04, 865.59) 


1.6 


(38.17, 748.19) 


(37.90, 742.79) 


(38.38, 752.39) 


(37.62, 737.19) 


(36.83, 721.43) 


(45.76, 899.99) 



TABLE II 

Time to synchronization [s](1st element of each 2-tuple) and energy consumption [10 _3 J](2nd element of each 2-tuple) under 

DIFFERENT PHASE RESPONSE FUNCTIONS (£j G (— 7T, 7r), I = gi — 92 — • - • — 9N — 0.01) 



eg \ £l 


0.05 


0.1 


0.2 


0.4 


0.8 


1.6 


0.4 


(22.93, 443.39) 


(23.17, 448.19) 


(23.14, 447.59) 


(22.58, 436.39) 


(21.53, 415.39) 


(22.44, 433.59) 


0.8 


(24.95, 483.79) 


(25.21, 488.99) 


(25.36, 491.99) 


(24.23, 469.39) 


(23.63, 457.39) 


(24.34, 471.59) 


1.6 


(30.14, 587.59) 


(31.92, 623.19) 


(31.75, 619.79) 


(30.35, 519.79) 


(28.15, 547.79) 


(29.09, 566.59) 



Using the same setup, we also implemented the network 
when the maximal/minimal £j is outside (— f , ■§)■ We ran 
the implementation for 100 times and each time chose the 
initial values of & randomly from a uniform distribution on 
(— 7r, tt). 28 of the 100 runs were unsynchronized, confirming 
Theorem [2] that all oscillators have to be connected to the 
global cue to guarantee synchronization. So we made g\ = 
92 = ■ ■ • = 9n = 0.01 and re-ran the implementation 
under different phase response functions. The results are given 
in Table II. With a decrease in e g , the energy consumption 
indeed decreases, confirming the effectiveness of our design 
methodology in Sec. [V] Moreover, using the same coupling 
strength, we also implemented the network under Peskin's 
phase response function used in [3], and obtained a (synchro- 
nization time, energy consumption) 2-tuple as (25.24, 489.59). 
Since it is larger than the smallest energy consumption in Table 
II, which is obtained under the same coupling strength, this 
confirms that by tuning the parameter e g in phase response 
function, energy consumption can indeed be reduced. 

Setting gi = g 2 = ... = .gjv = g, £ g = 0.4, and 
Si = 0.05, we also implemented the network under different 
strengths of global and local cues, i.e., different pairs of g 
and I. For each pair of g and I, we ran the implementation 
for 100 times, and each time we chose the initial values of 
£i randomly from a uniform distribution on (—tt, tt). The 
average time to synchronization is given by the first element 
of each 2-tuple in Table III. From Table III, we can see that 
a larger g indeed leads to a faster synchronization rate (a 
smaller synchronization time), whereas a larger I does not 
necessarily bring a faster synchronization rate. A larger I may 
even desynchronize the network when g is small (as illustrated 
by the last two elements of the first row). This confirms the 
analytical results in Remark [7] which state that the local cue 
may increase or decrease the synchronization rate. The same 
conclusion can be drawn for energy consumption, which is 
given by the second element of each 2-tuple in Table III. 

VII. Conclusions 

The synchronization rate of pulse-coupled oscillators is 
analyzed. It is proven that in addition to the strengths of 



global and local cues, the phase response function also de- 
termines the synchronization rate. This inspires us to increase 
the synchronization rate by choosing an appropriate phase 
response function when the phase response function is a 
design parameter. An application is the clock synchronization 
of wireless networks, to which pulse-coupled synchronization 
strategies have been successfully applied. By exploiting the 
freedom in the phase response function, we give a new design 
methodology for pulse-coupled synchronization of wireless 
networks. The new methodology can reduce energy consump- 
tion in clock synchronization. QualNet experiments are given 
to illustrate the analytical results. 

References 

[1] R. Mirollo and S. Strogatz. Synchronization of pulse-coupled biological 
oscillators. SI AM J. Appl. Math., 50:1645-1662, 1990. 

[2] A. Abdrabou and W. Zhuang. A position-based QoS routing scheme for 
UWB mobile ad hoc networks. IEEE J. Sel. Areas Commun., 24:850- 
855, 2006. 

[3] Y. W. Hong and A. Scaglione. A scalable synchronization protocol for 

large scale sensor networks and its applications. IEEE J. Sel. Areas 

Commun., 23:1085-1099, 2005. 
[4] A. Tyrrell, G. Auer, and C. Bettstetter. Emergent slot synchronization 

in wireless networks. IEEE. Trans. Mob. Comput., 9:719-732, 2010. 
[5] G. Werner-Allen, G. Tewari, A. Patel, M. Welsh, and R. Nagpal. Firefly 

inspired sensor network synchronicity with realistic radio effects. In 

Proc. SenSys 05, pages 142 -153, USA, 2005. 
[6] R. Pagliari, Y. W. P. Hong, and A. Scaglione. Bio-inspired algorithms 

for decentralized round-robin and proportional fair scheduling. IEEE J. 

Sel. Areas Commun., 28:564-575, 2010. 
[7] Y. Q. Wang, F. Nunez, and F. J. Doyle III. Energy-efficient pulse- 
coupled synchronization strategy design for wireless sensor networks 

through reduced idle listening. IEEE Trans. Signal Process., Accepted. 
[8] Y. Q. Wang and F. I. Doyle III. On influences of global and local 

cues on the rate of synchronization of oscillator networks. Automatica, 

47:1236-1242, 2011. 
[9] P. Delellis, M. di Bernardo, and M. Porfiri. Pinning control of complex 

networks via edge snapping. Chaos, 21:033119, 2011. 
[10] H. Kopetz and W. Ochsenreiter. Clock synchronization in distributed 

real-time systems. IEEE Trans. Comput., C-36:933-940, 1987. 
[11] E. Izhikevich. Dynamical systems in neuroscience: the geometry of 

excitability and bursting. MIT Press, London, 2007. 
[12] S. Barbarossa and G. Scutari. Bio-inspired sensor network design: 

Distributed decision through self-synchronization. IEEE Signal Process. 

Mag., 24:26-35, 2007. 
[13] F. C. Hoppensteadt and E. M. Izhikevich. Weakly connected neural 

networks. Springer, New York, 1997. 
[14] T. S. Rappaport. Wireless communications: principles and practice. 

Prentice Hall, New York, 2002. 



8 



TABLE III 

Time to synchronization [s] (1st element of each 2-tuple) and energy consumption [10 _3 J] (2nd element of each 2-tuple) under 

DIFFERENT STRENGTHS OF GLOBAL AND LOCAL CUES (£j £ (-IT, 7r), g\ = g 2 = ■ ■ ■ = 9JV = 9) 



9\l 


0.01 


0.02 


0.03 


0.04 


0.05 


0.06 


0.01 


(22.93, 443.39) 


(23.21, 448.99) 


(26.37, 512.19) 


(27.26, 529.99) 


(no sync, — ) 


(no sync, — ) 


0.02 


(17.49, 334.59) 


(18.90, 362.79) 


(22.03, 425.39) 


(24.60, 476.79) 


(24.38, 472.39) 


(21.39, 412.59) 


0.03 


(14.18, 268.39) 


(14.99, 284.59) 


(18.03, 345.39) 


(19.93, 383.39) 


(20.35, 391.79) 


(19.91, 382.99) 



[15] S. J. Park and R. Sivakumar. Load-sensitive transmission power control 
in wireless ad-hoc networks. In GLOBECOM, pages 42-A6, Taibei, 
2002. 

[16] C. V. Vreeswijk, L. F. Abbott, and G. B. Ermentrout. When inhibition 
not excitation synchronizes neural firing. J. Comput. Neurosci., 1:313- 
321, 1994. 

[17] H. K. Khalil. Nonlinear systems. Prentice Hall, New Jersey, 2002. 
[18] C. Godsil and G. Royle. Algebraic graph theory. Springer, New York, 
2001. 

[19] R. Horn and C. Johnson. Matrix analysis. Cambridge University Press, 
London, 1985. 

[20] P. Monzon and F. Paganini. Global considerations on the Kuramoto 
model of sinusoudally coupled oscillators. In Proc. 44th IEEE Conf. 



Decision Control, pages 3923 -3928, Spain, 2005. 

[21] S. Chung and J. Slotine. On synchronization of coupled Hopf-Kuramoto 
oscillators with phase delays. In Proc. 49th IEEE Conf. Decision 
Control, pages 3181-3187, USA, 2010. 

[22] A. Tyrrell, G. Auer, and C. Bettstetter. Fireflies as role models for 
synchronization in ad hoc networks. In Proc. Int. Conf. Bio Inspired 
Models of Network, Information and Computing Systems (BIONETICS), 
pages 1-7, Italy, 2006. 

[23] B. Ermentrout. Type I membrances, phase resetting curves, and 
synchrony. Neural Comput., 8:979-1001, 1996. 

[24] QualNet 4.5 User's Guide. Scalable networks inc. http://www. scalable- 
networks. com, 2008. 



